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The accuracy and physical significance of 
the classical Rayleigh-Sommerfeld and 
Kirchhoff difiraction integrals are assessed 
in the context of Sommerfeld's rigorous 
theory of half-plane diffraction and 
Maxwell's equations. It is shown that the 
Rayleigh-Sommerfeld integrals are in sat- 
isfactory agreement with Sommerfeld's 
theory in most of the positive near zone, 
except at sub-wavelength distances from 
the screen. On account of the bidirectional 
nature of diffraction by metallic screens 
the Rayleigh-Sommerfeld integrals them- 
selves cannot be used for irradiance calcu- 
lations, but must first be resolved into their 
forward and reverse components and it is 
foimd that Kirchhoff s integral is the 
appropriate measure of the forward irradi- 
ance. Because of the inadequate boundary 



conditions assumed in their derivation the 
Rayleigh-Sommerfeld and Kirchhoff inte- 
grals do not correctiy describe the flow of 
energy through the aperture. 
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1. Introduction 

In a previous paper [1] this author derived mathemat- 
ically rigorous expressions for the classical Rayleigh- 
Sommerfeld and Kirchhoff boundary- value diffraction 
integrals pertaining to circular apertures and slits illu- 
minated by normally incident plane waves. In spite of 
their functional differences, these diffraction integrals 
were found to be surprisingly similar and nearly indis- 
tinguishable in most of the near zone. They exhibited 
significant differences only in the immediate proximity 
of the aperture, but in this region their physical proper- 
ties were obscured by the fact that they or their normal 
derivatives, or both, do not reproduce the assumed inci- 
dent field. In these circumstances it was not possible to 
assess their physical significance by merely comparing 



them to one another In the present paper, they will be 
re-examined by applying them to the specific case of 
diffraction by a reflecting half plane and their physical 
properties will be interpreted in the context of 
Sommerfeld's [2] rigorous theory of half-plane diffrac- 
tion and Maxwell's equations. 



2. Comparison of Scalar Wave Functions 

The scalar wave functions U discussed in this paper 
all denote the complex disturbance at a point of ob- 
servation P(x, y, z) in the diffraction pattern of a per- 
fectly conducting, infinitesimally thin, semi-infinite 
screen that occupies the half plane x>Q, z = of a 
cartesian coordinate system, as depicted in Fig. Al of 
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Appendix A. The primary field is assumed to be a 
monochromatic plane wave with irradiance Eq, wave- 
length A, and circular wave number k = In/X that is 
normally incident from the half space z < and is 
plane polarized so that, in accordance with Maxwell's 
equations, dU/dz or U are continuous and equal to 
zero on crossing the screen. The resulting diffraction 
pattern is independent of y and will be denoted by 
U = yJEg u{x, z), so that we have |m| < 1 . 

Sommerfeld's half-plane theory dates back to the late 
1 800s and used to be discussed at length in textbooks 
[3-5]. However, it appears to be no longer included in 
modem curricula of theoretical optics, and therefore its 
main features are summarized and supplemented by 
new expressions for the diffracted irradiance in 
Appendix A, below. On combining Eqs. (A3a-c) and 
(A8a,b) of Appendix A it follows that, for normally 
incident light, Sommerfeld's solution is reduced to 



,(p.s) 



(X,Z)=Mg-FMg, 



:e*V(p), ;is=e*V(p), 



(la) 



(lb) 



p = J— [-sign(z)Vr + x-Vr-x], z 5^0, 

p = J— [-sign(z)Vr+x+Vr-x], Z5^0, (Ic) 



-4. 



where r = Vx + z and V(p) is the complex Fresnel- 
type integral defined by Eq. (A3d). These expressions 
are rigorously valid everywhere in the xz-plane of Fig. 
Al, except that along the x-axis p and p must be eval- 
uated as 



'Wt-^= 



P = -2sign(x), 



2sign(x)j^. r=-0. 



■2j^, z = 4^. 



(Id) 



(le) 



where z = ±0 refers to the positive and negative sides of 
the screen, respectively. This distinction is necessary 



because Wg andDwg /3z are discontinuous on cross- 
ing the screen, and is taken into account in 
Sommerfeld's theory by "wrapping" the diffracting half 
plane in a semi-infinite, two-sided Riemann surface so 
that its positive and negative sides are distinguished by 
the values 2k and of the polar angle ^ in Fig. Al . 

The corresponding results given by the Rayleigh- 
Sommerfeld theory are obtained from Eqs. (10a,b) of 
Ref. [1] by suitably modifying the limits of integration, 
leading to 



P = k4{^-xf+z\z>0, 



(2a) 



^RS V-^' ^) ~ 



ikz 



]d(^-x) 



nf^ip) 



, z>0, (2b) 



where Hj_^' = J^ +iY^ are Hankel functions of the first 
kind and Mth order. These expressions are valid for 
z > 0, only, and will be supplemented in this paper by 
the assumptions made in their derivation for z < 0; 
namely, 



a<'(jc,-fO) 



ik or 0, MRg(x, -tO) =1 or 0, 



and 



X < or X >0, for z = -tO, 



M^g'(x,z) = e or (e -i-e ), 
jc<Oor jc>0, for z<-0 



(2c) 



(2d) 



respectively. 

KirchhofFs diffraction integral, which will be required 
for the discussion in Sec. 3 is equal to the arithmetic 
mean of the Rayleigh-Sommerfeld integrals (2a,b), 

«K iX' z) = -[4s^ {x, z) + li^lix, z)], (2e) 

and can therefore be easily deduced from the above 
expressions. 

In the paraxial Fresnel approximation where z is posi- 
tive and large compared to X and xlz is small all of the 
above-mentioned solutions converge to the familiar 
Fresnel limit Mp. That is. 



,(p,s) 



Ul''' {X,z) ~ U^^'''{x,z) ~ Mj-(x,z) ~ Mp(x,z) 



-"^(-# 



z »j x\. 



(3a) 



where the right-hand expression follows from Eqs. 
(la-c) by letting ^r + x~4z [l±jc/(2z)], so that 

p = -x-Jll Xz, p = -yJSz/ X, and u^ « u^.^ The same 
result is obtained from Eqs. (2a,b) on replacing H*,^', 
Hf', and (3 by the leading terms of their asymptotic and 

For large negative values of z, Sommerfeld's theory also yields a 
comlementary term for the Fresnel approximation of the reflected 
diffraction pattern on the lit side of the half plane. 
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Taylor expansions. The Fresnel approximation, Eq. 
(3a), is estimated to be accurate within 1 % for z » 
lOOA. 

For numerical applications it is also useful to know 
that the above solutions all predict the same value, 



ut'\<d,z)=ut^'\<d,z)=u^{<d,z) 
= Mp(0, z) = — e ' , 



(3b) 



in the positive shadow boundary (x = 0, z > 0). In the 
case of the Rayleigh-Sommerfeld integrals, Eqs. (2a,b), 
this result follows from the identity 



Jd/H«(V7+7) = e'°, 



(3c) 



and was used in this work as the starting value for 
recursive numerical integrations as described in Ref 

[6]. 

The above expressions for Mg'''^ {x, z) and M^g'^' (x, z) 
were used to compute the squared magnitudes of these 
functions in the immediate proximity of the positive 
and negative sides of the aperture plane, as shown in 
Figs. 1 and 2. For these computations, Eqs. (2a,b) were 
evaluated as noted above and the Fresnel sine and 
cosine integrals required for the computation of V(p) 
and V(p) were evaluated using the algorithms of Ref 
[7]. The main conclusions drawn from these results are 
as follows. 

(1) On the positive side of the aperture plane the 
Sommerfeld and Rayleigh-Sommerfeld solutions are 
surprisingly similar, even at very small distances z. The 





l-PlAUiilHtfAl,^ - + ll.lii 



^X 



real and imaginary parts of Mg'''*' and M^g'^' contribut- distance z = +0.lA&om the aperture plane 



Fig. 1. I u'^''\x,z) I' ( ) and | u';''\x,z) \" ( ) vs x/X at the 



ing to the results plotted in Fig. 1 agree within + 1 % or 
better for z = 0.lA, and additional computations 
showed that this agreement improves rapidly for larger 
values of z. It follows that for all practical purposes the 
Rayleigh-Sommerfeld integrals are adequate for com- 
putations throughout the positive near zone, and hence 
it may be inferred that this will also be the case for the 
corresponding solutions for circular apertures and slits 
derived in Ref [1]. 

(2) The agreement for negative values of z is unsat- 
isfactory. In Sommerfeld's theory diffraction manifests 
itself as a field phenomenon that occurs on both sides 
of the aperture plane, so that the incident geometrical 
field is modified before it reaches the screen. On the 
other hand, in the Rayleigh-Sommerfeld theory diffrac- 
tion on the source side is explicitly ruled out, and here 
the results obtained from Sommerfeld's theory show 
that the assumed geometrical field (2d) is only a crude 



approximation of the true field. Thus, the main problem 
with the Rayleigh-Sommerfeld and Kirchhoflf integrals 
appears to be not so much that they fail to reproduce the 
assumed geometrical field values, but that the latter are 
themselves objectionable. 

The residual differences between u^/-' and u^''"'-' for 
z > can be attributed to the imperfect boundary condi- 
tions assumed in the Rayleigh-Sommerfeld theory. 
These boundary values are step functions that violate 
the wave equation and are the probable cause of the 
fact, shown in Appendix B, that the Rayleigh- 
Sommerfeld integrals also do not obey the wave equa- 
tion in the immediate proximity of the aperture plane. 
Although this wave-equation failure is small in most of 
the near zone, and thus unimportant for practical pur- 
poses, it is worthwhile to mention that it might be reme- 
died by replacing the boundary values Eq. (2c) with the 
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i«r 



^ — 

-5 



P'Paiarizaliim, t = - 0.1a 



p-Fokirization. z = + 

1 xlX 2 



\u\ 




s-Polanzalum , z = -0. 1 A 



xlX 



s-Poiarization, z = + 



xlX 



Fig. 2. I u-^-'\x,z) I' ( ) and | ul^\x,z) f ( ) vs xlX at the 

distance z = -0.11 from the aperture plane. 

corresponding values given by Sommerfeld's theory for 
z = +0; namely, 

l^""(-'^)=lor2V(-2,&-l,<Ae-'-!>, 

■\k dz ^ VA^ av2x 

X <0 or X >D, (4a) 

M<'^(jc,-K)) =1 or 2V(-2 /-), jc<0 or jc>0. (4b) 

The real and imaginary parts of these functions are 
plotted in Fig. 3, where it should be noted that du'f^ldz 
is discontinuous and singular, and Mg'' is not continu- 
ously differentiable, forx = 0. Nonetheless, they consti- 
tute improved boundary values because Sommerfeld's 
theory obeys the wave equation even at the diffracting 
edge itself (see Appendix B). 



Fig. 3. Real ( ) and imaginary ( ) parts of the boundary values 

Eqs. (4a,b) predicted by Sommerfeld's theory for p- and s-polarized 
lincident light. 



When Eqs. (4a,b) are substituted into the derivation 
of the Rayleigh-Sommerfeld integrals for the half plane 
one finds 

«<''>(x,z)^ Jd(^-x)^^^%i^H«(;8), z>0. 



dz 



(4c) 



«« (X, z) - ^ J d(,? - x)u'i\x, +0) ^^, z >0, 
2 ^ p 

(4d) 

where the integration now extends from -0== to +0°. 
Because the boundary values Eq. (2c) and Eqs. (4a,b) 
are the same for x < and the former are zero for x > 0, 
this can be rewritten as 
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,(p.^) 



(x, z) = M^g'*' (x, z) + Am^^ 



.") 



(4e) 



where 



Am 



(.p) _ „(p) 



RS 



Mq 



„(rt 



A i(fa-|) - 



= kj d(^ -x)[V(-2 ^) --^ jf e'""^' ]H«(i3), (4f) 

A<'=«<''-<^3'=i^^zJd(^-x)V(-2^)^^ 

(4g) 

are correction terms that can be added to the Rayleigh- 
Sommerfeld integrals to convert them to the exact val- 
ues given by Sommerfeld's theory. These expressions 
should be free of errors because Eqs. (4c) and (4d) are 
rigorous expressions of the Helmholtz' theorem in 
which Mg''^ and ulf-' are the same on both sides of the 
equal sign. 

This method was originally proposed by Braunbek 
[8-10], who envisioned its use for constructing 
improved solutions for large apertures of finite width 
and are bounded by straight or even curved edges. 
Braunbek's work involved the assumption that 
3Mg''V3z and Wg^' rapidly become negligibly small on 
the dark side of the screen, so that the effective ranges 
of integration in Eqs. (4f,g) are only a few wavelengths 
wide and approximative methods can be used. 
According to Fig. 3 this is a valid assumption for 
3Mg''V3z but not for u^'\ so that computational difficul- 
ties could be encountered in the case of\M^g . 



ponents in the forward and reverse directions. For nor- 
mally incident light, these respective directions are par- 
allel and anti-parallel to the unit vector n = [0,0,1] in 
the direction of the positive z-axis, and the final expres- 
sion for the time-averaged Poynting vector (A7c) is 



<Ss> 



[EJx,z)-EJx,z)]n = E. 



Mg |')«, (5) 



where E^ and E^ are the forward and reverse irradiances 
incident on the opposite sides of any given area element 
dx dy containing the point of observation P.'' These irra- 
diances are given by the squared magnitudes of the 
basic Sommerfeld fiinctions Mg and Mg themselves, and 
thus the quantities |Mg - Mgp or |Mg + Mgp do not represent 
the irradiances of the field forp- and s'-polarized light. 
Accordingly, the forward and reverse irradiances of the 
field are independent of the state of polarization of the 
incident light, and in this connection it should also be 
noted that in practice the reverse irradiance E^ is not 
easily observable as it may be obscured by a detector 
placed in the path of the forward field. 

It now seems reasonable to interpret the Rayleigh- 
Sommerfeld theory in a like manner, so that the quanti- 



ties M^g and M^g defined by Eqs. (2a,b) are also regard- 
ed as bidirectional wave fiinctions that can be resolved 
into mutually incoherent forward and reverse compo- 
nents, Mr and Mr. Thus we define, in analogy to Eq. (la). 



"k = -[Hrs(^' z) + M^s^(x, z)], 



Mj, L'^RS (^J-^) '^RSV^'-^/J- 



(6a) 



3. Irradiance and Energy Flow 

Although the squared magnitudes of scalar wave 
fiinctions are commonly identified with the irradiance 
of the field, the data plotted in Figs. 1 and 2 must not 
be interpreted in this manner The diffracted field spec- 
ified by Sommerfeld's solution is a bidirectional field 
composed of two plane waves, Mg and +Ms which prop- 
agate in the opposite directions of the incident primary 
field and its reflection from the screen. When 
Maxwell's equations are invoked, as in Eqs. (A5) 
through (A7) of Appendix A, it is found that in accor- 
dance with the principle of interference these waves 
cannot interfere with one another^ so that the effective 
energy flow is composed of mutually incoherent com- 



Except for grazing incidence, where interference can occur as in 
Fresnel's mirror experiment. 



and hence it follows that the corresponding forward and 
reverse irradiances, E^^ and E^ will be given by an 
expression analogous to Eq. (5), 



<S^>- 



[E^ix,z)-E^(x,z)]n = E^(\ u^f -\u^ |')«. 

(6b) 



It will be noted that the forward wave function Uf^ 
defined by Eq. (6a) and Kirchhoffs integral (2e) are 
identically the same, and therefore the subscript "K" 
was retained in the above equations. The Kirchhoff and 
Rayleigh-Sommerfeld solutions were originally 
derived on the mutually exclusive assumptions of black 



For oblique incidence of the primary field, as assumed in Eq. (A8a), 
these irradiances are defmed with respect to area elements which are 
perpendicular to the directions of propagation of the incident and 
reflected geometrical field. 
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and metallic screens, and it is generally agreed that Eq. 
(2c) has no definable physical meaning as it would 
somehow imply the coherent superposition of two 
orthogonal states of polarization. However, in the pres- 
ent context, the Rayleigh-Sommerfeld integrals are 
interpreted as composite quantities and taking their sum 
and difference is tantamount to resolving them into 
their basic components. Accordingly, Kirchhoffs inte- 
gral Mf^now appears as an integral part of the Rayleigh- 
Sommerfeld theory for metallic screens so that m*^^ and 
M^5 provide the framework for the evaluation of all 
field parameters while m^ and its counterpart m^ define 
the flow of field energy. This new interpretation of 
Kirchhoffs integral has a precise, physically realizable 
meaning. 

A numerical comparison of the forward irradiances 
E^ and E^^ defined by Eqs. (5) and (6b) is presented in 
Figs. 4 and 5. As expected, these quantities are essen- 
tially the same on the positive side of the aperture 
plane, the agreement being on the order of a few per- 
cent for z = +0. 1 A and increasingly better for larger val- 
ues of z. This confirms that the identification of \u-gf 
with the forward irradiance £k is a valid assumption. As 
also expected, the agreement is poor on the negative 
side because in this region £,. represents only the undif- 
fracted geometrical field. The even symmetry of the 
irradiance E^ shown in Fig. 5 suggests that the modifi- 
cation of the geometrical field due to diffraction is 
isotropic in the immediate vicinity of the edge. 



4. Conclusions 

The above comparison of the classical Rayleigh- 
Sommerfeld boundary-value theories with 
Sommerfeld's rigorous theory for diffraction by a per- 



aaaA 



.0 I 




Fonvard Irradiance. z^-0 \X 



xlX 



xlX 



Fig. 4. Forward irradiances E^{x^) ( ) and -Bk(-'^,z) ( ) vs x/1 at 

the distance z = +0.11 from the aperture plane. 



Fig. 5. Forward irradiances E^{x,z) ( ) and -Bk(-'^,z) ( ) vs xlX at 

the distance z = -0. 11 from the aperture plane. 

fectly reflecting half plane has added substantially to 
the understanding of the physical significance of these 
theories. 

It was found that the mathematical expressions and 
algorithms presented in Ref [1] for the Rayleigh- 
Sommerfeld integrals are in very satisfactory agree- 
ment with Sommerfeld's half-plane theory. Thus, they 
are well suited for computations in most of the positive 
near zone, and it is inferred that this will also be the 
case for the corresponding Rayleigh-Sommerfeld inte- 
grals and slits derived in Ref [1]. Sommerfeld's theory 
also confirms that, on the whole, the differences 
between these respective solutions for p- and s'-polar- 
ized incident light are small so that polarization effects 
are small, as might be expected for normally incident 
light. All in all, it appears that the use of Helmholtz' 
theorem has proved remarkably effective in compensat- 
ing for the inadequate boundary conditions assumed in 
deriving the classical boundary-value integrals. The 
residual differences between the Rayleigh-Sommerfeld 
and Sommerfeld solutions are confined to sub-wave- 
length differences from the screen, and it is shown in 
Appendix B that in this region the former do not obey 
the wave equation. 

The comparison with Sommerfeld's theory and its 
interpretation in terms of Maxwell's equations has also 
revealed a previously overlooked aspect of diffraction 
by a reflecting screen; namely, that the optical field is 
bidirectional and comprises light traveling in opposite 
directions even on the positive side of the screen. 
According to the principle of interference, the observ- 
able Poynting vector is given by the incoherent vector 
sum of its components in the forward and reverse com- 
ponents, and thus it is impermissible to express the 
near-zone irradiance of the field as the squared magni- 
tudes of scalar wave functions. Rather, the latter must 
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be resolved into their forward and reverse component 
and it turns out that Kirchhoff s integral is the appropri- 
ate expression for the forward irradiance of the field 
even in the Rayleigh-Sommerfeld theory. The forward 
and reverse irradiances were found to be independent 
of the state of polarization of the incident field. 

It was noted that the residual deficiencies of the 
Rayleigh-Sommerfeld and Kirchhoff solutions in the 
proximity of the positive aperture plane can be 
removed by replacing the originally assumed boundary 
values with those predicted by Sommerfeld's theory. 
This was not be pursued further as it would produce 
only marginal improvements on the positive side of the 
screen, without removing the problem that the classical 
boundary-value integrals all exhibit discontinuities 
with respect to the incident geometrical field. A more 
effective approach would be the derivation of improved 
approximations for the entire field by constructing ana- 
lytical continuations of the existing boundary-value 
solutions into the half space z < 0. This will be attempt- 
ed in a subsequent publication. 



5. Appendix A. 
Theory 



Sommerfeld's Half-Plane 



(|) = 71 




([) = a 



([) = 7i+a 



P(r,(^) 



^ = 2k 



^ = 



X 




(|) = 7r™a 



Fig. Al. Basic geometry and notation for Sommerfeld's theory. 



In Sommerfeld's rigorous treatment of diffraction by 
a straight edge the screen is assumed to be a perfectly 
conducting, infinitesimally thin, semi-infinite sheet that 
covers the half-plane x > 0, z = of the Cartesian coor- 
dinate system in Fig. Al . It is assumed, further, that the 
primary field is a monochromatic plane wave, p- or s- 
polarized with respect to the xz-plane and incident upon 
the screen in a given angular direction a. As the optical 
field so defined must be independent ofy it will be con- 
venient to use cylindrical coordinates (r, (/>,>>) given by 



=V^ 



r = ylx + z , (j) = arccos(— ), x = rcos (j), z = -rsin (j), 
r 

(Ala) 

so that (/> and a are measured clockwise from the posi- 
tive X-axis, and the illuminated and shaded sides of the 
screen are distinguished by (/> = and 2%, respectively.'* 
In this notation, the primary field, its reflection by the 
screen, and the unit vectors in their respective direc- 
tions of propagation are 



For numerical computations the expressions = ACOS(x/r) if z < 
and = 27t-ACOS(x/r) if z < must be used, where ACOS(x/r) is the 
principal value of arccos {xir). 



U,^.=^E,c-'^^^''"\ ^e„. =+V^e-'*'-^<*-'. (Alb) 

where the time factor of the field is assumed as e""*, 
k = 2tiIX is the circular wavenumber of the light, £o is 
the incident irradiance, and the dual sign of ^geom 
accounts for polarization-dependent phase changes on 
reflection. 

According to these definitions, the diffracted field at 
a given point P(r,(/>) must obey the scalar wave equa- 
tion, 



^U 



(p.s) 



^^yip.") 1 ^i/ip.") 1 ^^iji 



I^TliP,") 



dr' 



r dr r^ 3(/>^ 



(A2a) 



as well as the boundary conditions for (/> = and (j) = 2%, 



t/*"' = or 



3z 



= 0, for 0=0 and 2 n:, (A2b) 



according as the light is p- or 5-polarized. Furthermore, 
in the limit r ^ oo these solutions must correspond to 
the optical field according to geometrical optics, and in 
this respect it is necessary to distinguish three regions 
of space as indicated in Fig. Al: 



63 



Volume 108, Number 1, January-Febraary 2003 

Journal of Research of the National Institute of Standards and Technology 



I. The reflection space (0 < < a), where the inci- 
dent and reflected waves are both present and the geo- 
metrical field is f/<f^' = f/^_ + ^_ . 

II. The transmission space (a < < a + t:), where 
only the incident wave is present and the geometrical 
field is [/<"■'' = U . 

geom geom 

III. The shadow space (a + t: < < 2%), where the 
geometrical field U^^^^ is zero. 

Sommerfeld's solution of the diffraction problem so 
defined is: 



Ul''-'\r,^) = 4E~,(u, + {^), 



(A3a) 



M„ = e 



-iArcos(0-a)T 



V(p), M5 =e-'*"°^**+«'V(p), (A3b) 



P = ^|cos 



1 



-(0-a) 



'"\'r"^ 



i+a) 



(A3c) 



e'" "f . ^ h 



V(p)=^JdTe^ =^l=C(p)+S(p)] 



-^C(p)-S(p)], 



(A3d) 



where C(p) and S(p) denote the usual Fresnel cosine 
and sine integrals. 

Although the derivation of these expressions is too 
complicated to be included in this paper, it is not diffi- 
cult to verify that they have the following properties: 

(1) As shown in Appendix B, Wg and Kg obey the wave 
equation, Eq. (A2a), everywhere in space, inclusive of 
the diffracting edge itself They represent plane waves 
which propagate with space-dependent amplitudes, 
V(p) and V(/5), in the respective directions of the unit 
vectors, Eq. (Ale). Except in the reflection space and at 
small distances r from the diffracting edge. Kg is signif- 
icantly smaller than u^, and in the limit r — > 00 Eqs. 
(A3a,b) are reduced to the above-mentioned geometri- 
cal solutions. 

(2) In addition to the usual plane-wave components 
which are proportional to Mg and Mg themselves, the 
derivatives, 



3Mg / dx 
3i5g / dx 



'3Mg /dz 
3i5g / dz 



= -iA:cosa 



+ COS 



r, 



:+iA:sina 



-sm 



^(0±«) 



-(0±«) 



,'<-!' 



Mr 
(A4a) 



;<-!' 



Mr 
(A4b) 



contain terms which involve the derivatives of V(p) 
and V(p) with respect to r and, thus, are singular as 
1/vr at the diffracting edge. This suggests the exis- 
tence of cylindrical waves which originate at the edge, 
but nonetheless the edge does not radiate energy. 
Although the radiant intensity /of the edge is infinite as 
1 / vr , the radiant flux MQ emitted into any given solid 
angle element dQ = rd0 dy is zero as r — > 0. These 
cylindrical waves are evanescent and vanish in the limit 
r»A. 

(3) Equations (A3b) and (A4a,b) show that (Kg - Kg) 
and 3(Mg - u^ydx are continuous on crossing the screen 
(0 = — > = 27:), whereas 3(Mg - Mg)/3z is not. 
Conversely, 3(Mg + u^ldz is continuous, whereas (Mg - 
Kg) and 3(Mg + u^ldx are discontinuous. These are the 
expected properties of the tangential (magnetic or elec- 
tric) field vector components when the light field is 
polarized parallel or perpendicular to the diffracting 
edge, and thus it is permissible to apply Maxwell's 
equations in the form 

H'"' = J2£^4-(M3-Mg)m, E'"^ =-,pcurlH<'", 

(A5a) 



-Hi 



E« = ^2f„^^(H,+^,)m.rf 



' --- '^curU^"', 



(A5b) 



where m = [0,1,0] is the unit vector in the direction of 
the positive 3^-axis, £ and // are the dielectric constant 
and magnetic permeability of the medium of propaga- 
tion, and normalization factors were used so that the 
squared field vectors have the dimension of irradiance 
[W/m^]. On substitution of the derivatives Eqs. (A4a,b) 
and re-arranging terms, this leads to 

E*''' = /2£ *\— \ -sma(« +u ), 0, -cosa(M -u ) 



0, cos 



a ] I (p cc 



fa] I f cc 

— H — -cos 

2 2 2 2 



2£„ 



uA + 



iky/Xr 



Ht + 



ik^Xr 



(A6a) 
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H^^' =J2Eil — <-sma(u +u ), 0, -coia{u -u ) 
^ "y^i L s s s s J 



\k\lXr 



a \ .10 a 

sin| — H — +sm 

2 2 I 2 2 



1 T « I ((pa 

0, cos — H — + cos 

2 2 2 2 



«^-t* 



uA- 



i^t- 



ikyJXr 



(A6b) 



ikyJXr 
where 

t = [-sina,0,-cosa], t = [sina,0, -cosa] (A6c) 
are the unit tangent vectors indicated in Fig. Al, and 



. , (j) a I _ ((j) a 
sm|- + - O,cos-+- 



. .(j) a] ((j) a 

sin ,0,cos 

'22 22 



(A6d) 



are unit vectors in the directions of the above-men- 
tioned evanescent cylindrical waves. The middle por- 
tions of Eqs. (A6a) and (A6b) are equivalent to the 
expressions cited by Bouwkamp [11] for the electro- 
magnetic field components given by Sommerfeld's 
solution. However, it appears that their representation 
in terms of the unit vectors defined by Eqs. (A6c,d) has 
not previously appeared in the literature. 

(4) Equations (A5a,b) and (A6a,b) express each of 
the electromagnetic field vectors as the sum or differ- 
ence of two components which, like Sommerfeld's 
wave functions u and u themselves, are easily recog- 
nized as representative of a forward or reverse wave 
motion. It is obvious that these components must be 
mutually incoherent because, in all cases of practical 
interest, they propagate in opposite or nearly opposite 
directions, thus precluding any interference between 
them. Therefore, the corresponding Poynting vectors 
are given by the expressions 



(*--!) 



E^"^ X H^"^* = 2£'„ { [-Ms t + ^-= v] 



[k\IXr 



i(i.4) 



* « e * ^ 

xu^m]-[u^t-i r^vJXMgin] |, 



(A7a) 
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£(») X H<'>* = 2E^ {ms m X [m* t - ^ v] 



+ Mgmx[Mgt- 



-i(*.-f) 



[k\ Xr 



V]}. 



(A7b) 



where it must be also be taken into account that, as 
mentioned above, the terms in l/^Jr have no energetic 
significance. Thus, these terms must be ignored and the 
following simple result is obtained for the observable 
Poynting vectors [5] of the field, 



<S>=-Re[E<'"xH<'"*]=-Re[E<^'xH<^'*] 
2 2 



(A7c) 



where E^ and E^ are the forward and reverse irradi- 
ances^ incident on area elements normal to the unit vec- 
tors in the directions of propagation of the incident and 
reflected field. 



s = mxt=[-cosa,0,sina], 
s = mxt=[-cosa,0, -sina], 



(A7d) 



as indicated in Fig. Al. These equations are not a part 
of Sommerfeld's original theory and may be interpreted 
as follows: 

(a) In the vicinity of the diffracting edge the wave 
functions Mg and Mg , and thus the forward and reverse 
irradiances E^ and E^ as well, are similar in magnitude. 
Therefore, light traveling in both directions is present 
on both sides of the screen except, in the Fresnel 
approximation where Kg is negligibly small. This is a 
significant departure from the classical formulation of 
Huygens' principle, where a reverse flow of energy on 
the positive side of the aperture is precluded by the 
explicit assumption that light does not travel back- 
wards. Instead, Sommerfeld's theory asserts the pres- 
ence of a bidirectional flow of energy on both sides of 
the aperture plane, and in this connection it is relevant 
to cite two earlier papers on Sommerfeld's theory by 
Braunbek [12] and Braunbek and Laukien [13]. The lat- 
ter includes an interesting diagram depicting a swirling, 
bidirectional flow of energy at sub-wavelength dis- 



In order to conform to established nomenclature, the letter "E" is 
used in this paper to denote irradiance as well as electric-field vec- 
tors. There should be no confusion, as these different quantities are 
consistently indicated in lightface italic and boldface roman type, 
respectively 
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tances from the diffracting edge. Although it appears 
that Braunbek and Laukien assumed a coherent super- 
position of the forward and reverse fields, the eddy cur- 
rents shown in their diagram can be regarded as 
Maxwellian analogues of Huygens' wavelets. 

(b) In spite of the explicit assumption of separate 
boundary conditions for p- and s'-polarized incident 
light, the forward and reverse irradiances E^ and E^ 
defined by Eq. (A7c) are the same in both cases. The 
composite wave functions u^ + u^ pertain to different 
states of polarizations only insofar as their phases are 
concerned, but their squared magnitudes cannot be 
used to describe the energy flow in the field as they 
contain non-observable cross terms in u^u^. This dis- 
tinction disappears in the Fresnel region, where the 
field on the positive side of the aperture plane is unidi- 
rectional (ms « Ms) and the usual definition of irradi- 
ance as the squared magnitude of the total wave func- 
tion is justified. 

(6) For a normally incident field (a = 7t/2) one finds 



cos((/> + a) =+sin0, 



(A8a) 



cos 



-(0-i-a) =-^ cos — ±sm — 

_2 J V2l^ 2 2 

- — U\+cos(p ±^l -coscj) ) 



(A8b) 



6.1 Sommerfeld's Solution 

In this subsection, Sommerfeld's solution (A3a-d) 
will be written as 



Ms orMs=M=e*"'''V(p), 

Porp = JyCos -^ I ^=(i> + a. 



Thus, 3/30 = 3/3)8 and 



3V 1 Tl r>^ i[fa-(l-cos/J)-|] 

^— = -^=cos -;8 e ■* , 

3V 1 ( ^ n\ i[*'(l-cos/J)-|] 
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30 4h- U ' 
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2r 2 



cos I — 8 
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(A9a) 



(A9b) 



— = -i^cospM+-;^cos -p b ■* , (A9c) 
3r 4Xr yi 



(A9d) 



Hence, the starting equations in Sees. 2 and 3 of the 
main text are obtained by using Eqs. (Ala) to reintro- 
duce Cartesian coordinates and noting that the sine 
term in the central portion of Eq. (A8b) is always posi- 
tive, while the cosine term and z are opposite in sign. 



6. Appendix B. 
Conformance 



Wave-Equation 



It is commonly agreed that one of the most important 
measures of the physical significance and mathematical 
rigor of scalar diffraction theories is whether, or how 
well, they satisfy the wave equation [Eq. (A2a)]. In this 
appendix, this aspect of the Sommerfeld and Rayleigh- 
Sommerfeld theories is analyzed. It is shown that 
Sommerfeld's solution obeys the wave equation rigor- 
ously everywhere in space, whereas the Rayleigh- 
Sommerfeld integrals (and, thus, Kirchhoff s integral as 
well) exhibit deviations from the wave equation in the 
immediate proximity of the aperture plane. 



l---^"-^-fl^V-^' 



3 u 



-{k^r^ sin^ P -ifo-cos j8)m 



'' ^•,2-2(1 

2 



cos 



a^)' 



(A9e) 



(*--) (A9f) 



Hence it follows immediately that the wave equation, 
Eq. (A2a), is rigorously satisfied. It should be noted 
that this is true everywhere in space, even at the dif- 
fracting edge itself where each component of the 
Laplace operator is singular. 

6.2 Rayleigh-Sommerfeld Integrals 

In this next subsection, the wave equation 
conformance of the Rayleigh-Sommerfeld integrals 
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Mrs"^ (x, z) is analyzed by computing numerical values 
of the quantities, 



</>(x,z) = 



1 






+1 



il^'\x,z), (AlOa) 



w^^^{x,z): 



M«(i3j 



]d(^-x) 



^vm 



H«(i3) , ikzfnfiP) 



(AlOh) 



which would everywhere be identically equal to zero if 
the wave equation is satisfied. For this purpose, the first 
and second derivatives of the Hankel functions appear- 
ing in Eqs. (2a,b) are evaluated by means of the identi- 
ties 






_d_ 



H«(i3) 



Hfm 



(A 10b) 
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di3 
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Hfip) dp__k^ dp _ k'z 
P' ' dx~ p ' dz~ p' 

(A 10c) 



Hence, the required jf-derivatives of m^^'^' are found by 
differentiating Eqs. (2a,b) with respect to the upper 
limit of integration, yielding 



a<'_m«(i3j du^^l jk'zRfHPJ 



dx 



2 dx 

P^=k^x^i 
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(AlOd) 



a^< _ k'xuf^pj 3X1 _ ^k'xzufipj 
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(AlOe) 
and the z-derivatives are obtained by substitution of 



3X'(ig) _ J^2 
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into Eqs. (2a,b). The final expressions are 



<S'(X,Z): 



ik'xzK'i\pJ 



ikz 



]d(^-x) 



H«(i3) 3H«(i3) , (kzfnf^p) 
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The integrals on the right-hand side of these expres- 
sions were evaluated by recursive numerical integration 
as described in Ref [7]. In accordance with Eq. (3b), 
the starting values used for these computations were 
w<r'(0,z) = 0. 

The real and imaginary parts of (A10h,i) are plotted 
in Fig. A2 for z = 0. 1 , showing that at this distance the 




MX 




Fig. A2. Real ( ) and imaginary ( ) parts of the quantities 

defined by Eqs. (A101i,i) at the distance z = +0.11 fiom the aperture 
plane. 
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wave-equation failure is substantial. Additional compu- 
tations indicated that the corresponding values of 
w'^"^ (x, z) decrease at larger distances z, but fall below 
0.01 only when z>30X so that calculable departures 
from the wave equation are present throughout the near 
zone. This is not surprising because, otherwise, the 
Rayleigh-Sommerfeld and Kirchhoff boundary-value 
integrals would be rigorously correct. On the other 
hand, the above results are at odds with a fallacious 
belief that these integrals must obey the wave equation 
because their integrands do. The fact of the matter is 
that the assumed boundary conditions Eqs. (2c,d) 
abruptly truncate the incident field at the edge and, 
thus, transform these integrands into discontinuous 
functions that violate the wave equation. 
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